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Abstract 

The production in hadron-hadron collisions of jet pairs with large rapidity separa- 
tion and comparable modest transverse momentum is, in principle, described by the 
perturbative QCD BFKL equation. The measurement of such jet pairs appears well 
suited to the LHC detectors, with their ability to detect forward jets. We present pre- 
dictions for dijet cross sections and correlations obtained using a BFKL Monte Carlo 
which allows kinematic and other subleading effects to be incorporated. The enhanced 
phase space for gluon emission at the LHC makes the BFKL behavior somewhat easier 
to observe than at the Tevatron, although kinematic effects are still important. The 
production of forward jets in association with heavy Higgs bosons via the gauge boson 
fusion mechanism is also studied and compared with QCD dijet production. 



1 Introduction 



In perturbative QCD physical quantities are typically computed in fixed-order expansions in 
powers of the coupling constant as- In some kinematic regimes, however, large logarithms 
multiply the coupling and the series must be resummed to all orders. In certain cases this 
resummation is possible and meaningful predictions can still be obtained. For example, large 
logarithms in the small-x region in deep inelastic lepton-hadron scattering are resummed by 
the Balitsky, Fadin, Kuraev and Lipatov (BFKL) equation |l|], to yield the well-known BFKL 
prediction F2 ~ (l/x)^, with A = 4Ca In 2 0:5/71 ^ 0.5. 

In hadron collisions BFKL resummation also applies to production of jet pairs at large 
rapidity separation A when the jets have comparable transverse momenta that are small 
compared to the center-of-mass energy. At large values of A, the cross section receives con- 
tributions from real and virtual gluons emitted in the rapidity interval between the two jets, 
as illustrated in Fig. ^(a). The BFKL equation resums the leading logarithmic contributions 
from these gluons 0, giving rise to a cross section that increases with A, d ~ exp(AA), 
with A defined above. Mueller and Navelet pointed out the potential of this increase, 
which does not appear in lowest-order QCD, to provide a way to investigate BFKL physics 
at hadron colliders by looking at dijet production at large rapidity separation. 

Unfortunately the A increase disappears when the subprocess cross section is folded in 
with parton distribution functions (pdfs), which decrease with A more rapidly than a in- 
creases. Several methods have been suggested to avoid this pdf sensitivity when studying 
'Mueller-Navelet' jets. One [0, ||, ^ is to study the decorrelation in the two jets' azimuthal 
angles that arises from emission of gluons between the jets; BFKL predicts a stronger decor- 
relation than does fixed-order QCD, and this prediction should be relatively insensitive to 
the pdfs. Another ^ is to look for the A increase in the ratio of dijet cross sections at 
two center-of-mass energies for corresponding parton momentum fractions, so that the pdfs 
cancel. 

In practice it is not useful to compare analytic asymptotic BFKL predictions directly 
with experiment because nonleading corrections can be large. In particular, in the analytic 
BFKL calculation gluons can be emitted arbitrarily, with no kinematic cost, and energy and 
momentum are not conserved. In Refs. ^ (see also [0]) a Monte Carlo approach is used in 
which the emitted gluons are subject to kinematic constraints, and other nonleading effects 
such as the running of as are included. Kinematic constraints are seen to have a significant 
effect, suppressing the emission of large numbers of energetic gluons. The studies in and 
0] focused on the azimuthal decorrelation and cross section ratios at the Fermilab Tevatron 
pp collider. 

Here we are concerned with dijet production at the CERN LHC pp collider. The phase 
space available for studying dijet production is larger at the LHC than at the Tevatron, 
because the center-of-mass energy is much higher while measurable jet transverse momentum 
thresholds increase only slightly, if at all. As a result we expect kinematic constraints to have 
weaker effects at the LHC, so that the predictions of naive BFKL should be more robust 
than at the Tevatron. We will address this in more detail below. 

BFKL physics (or, more generically, QCD) is not the only potential source of jet pairs 
with large rapidity separation at the LHC: heavy Higgs production via gauge boson fusion in 
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the process qq qqH (see Fig. |I|(b)) is well known to give rise to forward jets 0. In fact 
the forward calorimetry in the LHC detectors has been designed to be able to detect jets at 
large rapidity for the purpose of observing this process |10|. Detection of jets with rapidity 
\y\ < 5 (i.e. A < 10) and transverse momentum px > O{20 GeV) should be possible. 

In contrast to the BFKL case, where the region between the jets is populated by gluons, 
the jets produced in association with the Higgs are characterized by an absence of hadronic 
activity between them — a 'rapidity gap.' We will see other differences when we compare 
other features of dijet production in the two cases below. The Higgs process therefore acts 
useful comparator for the BFKL study. 

The paper is organized as follows. In the next section we summarize the essential features 
of the BFKL MC calculation used to obtain our results. In section 3 we present results 
for dijet production in pp collisions at center-of-mass energy 14 TeV, first in the BFKL 
approach, and then for jets produced in association with a heavy Higgs boson. In section 4 
we investigate the collision energy dependence of the BFKL dijet cross section and present 
the ratio of cross sections at 10 and 14 TeV. In section 5 we present our conclusions. 



2 Improved BFKL approach to dijet production in had- 
ron collisions 

The differential subprocess cross section for production of a jet pair with rapidity difference 
A = — 1/2I and relatively small transverse momenta Pti,Pt2 is given in the BFKL approach 
by 

— "^^^^ " A) (1) 

The Laplace transform of the function / with respect to A satisfies the BFKL equation. 
This equation contains integrals over both real and virtual gluons, and when its Laplace- 
transformed solution is substituted into Eq. (|l|), it gives the sum over emitted gluons and 
total cross section behavior discussed above. This is the 'naive' BFKL case. 

As indicated in the introduction, the price of having an analytic solution is integration 
over radiated gluons with arbitrarily large transverse energies. We can improve on this 
approach by unfolding the implicit sum over gluons that appears in the analytic solution. 
Details can be found in ||^; see also for a similar approach. The basic idea is that real 
gluons are separated into 'resolved' and 'unresolved' gluons, the latter being those that fall 
below a transverse momentum threshold /xq and are undetectable as jets. The unresolved 
gluon contributions are then combined with those from virtual contributions to give an 
overall form factor that suppresses the cross section. The cross section is then an explicit 
sum over resolved real gluons; for fixed as it becomes 
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with 7^'^*') = 5{pti + PT2) and T^'-"^ for > given in [^], where the expressions for running 
as can also be found. The analytic result for the total cross section is reproduced in the 
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limit /io — > 0.0 

It is straightforward to implement this in a Monte Carlo program where kinematic cuts 
can be applied directly. The total cross section is obtained by folding in the subprocess cross 
section with the appropriate parton distributions functions; see for details. 



3 Jet pair production in pp collisions at ^/s = 14 TeV 
3.1 BFKL calculation of dijet production 

The results shown here were obtained using the BFKL Monte Carlo program described above 
and in Ref. ^j. We consider dijet production in which the jets have equal and opposite 
rapidity yi = —1/2 = A/2 and transverse momentum greater than some value -Pr-Q For the 
LHC we consider transverse momentum thresholds Pt = 20 GeV and Pt = 50 GeV. Given 
the high level of hadronic activity at the LHC, the former value may be slightly optimistic 



but the latter is certainly achievable [10|. We use CTEQ4L parton distributions ||1 1|| . 

Figure ^ shows the dijet production cross section at the LHC with center-of-mass energy 
14 TeV for pt > 20 GeV (upper curves) and pt > 50 GeV (lower curves). All of the curves 
fall off with increasing A as a result of pdf suppression; note that the difference in momentum 
thresholds causes the pt > 50 GeV cross sections to reach the kinematic limit at smaller 
values of A than the px > 20 GeV curves. 

The naive BFKL prediction, obtained using the analytic (asymptotic) BFKL solution 
combined with pdfs assuming lowest-order kinematics, is shown as dashed curves in the 
figure. For both pjn thresholds, the naive BFKL cross section is largest for all values of A; 
this reflects the rise with A of the subprocess cross section. The falloff with A is due to the 
pdfs. The BFKL MC results, shown as solid curves, lie well below naive BFKL as a result of 
the kinematic suppression.^ Clearly nonleading effects due to kinematics are important. Also 
shown for comparison are the cross sections in the asymptotic (large A) lowest-order QCD 
limit (see for details about this limit). It lies below the naive BFKL curves, because the 
subprocess cross section is asymptotically fiat in A, and it lies above the BFKL MC curves 
because it reflects only leading-order kinematics, whereas the BFKL MC cross section pays 
the pdf penalty for emitted gluons. 

A characteristic of BFKL dijets that clearly distinguishes them from those in leading- 
order QCD is a decorrelation in azimuthal angle at large values of A 0, ||, ||. In QCD 
at lowest order the two jets are strictly back-to-back in the transverse plane, so that their 
azimuthal angles are completely correlated, independent of rapidity. In contrast, in the 
BFKL case gluons are radiated in the rapidity interval between the two jets, and as this 
interval increases more gluons can be radiated. The presence of these gluons weakens the az- 
imuthal correlation between the two jets of interest, so that one sees an increasing azimuthal 
decorrelation as A increases. 



^In practice, the value ^0 = 1 GeV gives a very good approximation to the hmiting result and acceptable 
Monte Carlo efficiency. In general /iq should be small compared to pti,PT2- 

^Equal and opposite rapidities are taken for convenience only and are not required for validity of the 
BFKL approach; the rapidity separation is the important quantity. 

■^The running of as also contributes to the suppression but kinematics is the main effect. 
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As discussed in [Q, applying kinematic constraints as in the improved BFKL MC gives 
less of a decorrelation than predicted in naive BFKL because emission of arbitrarily many 
gluons of arbitrary energy is no longer allowed. In fact, for any given transverse momentum 
threshold, there is some A at which the jet pair alone saturates the kinematic limit, and 
emission of additional (real) gluons is completely suppressed and the azimuthal correlation 
returns. As we approach that limiting value of A we therefore expect to see a transition 
back towards correlated jets. 

This is illustrated in Fig. ^, where we show the azimuthal decorrelation in dijet production 
in the improved BFKL MC approach (upper curves). We show the mean value of cosA0 
where we have defined 

A0= 101 -02| -TT (3) 

SO that A0 = when the two jets are back-to-back in the transverse plane. The jets 
are completely correlated at A = 0, and as A increases we see the characteristic BFKL 
decorrelation, followed by a flattening out and then an increase in (cos A0) as we approach 
the kinematic limit. Not surprisingly, the kinematic constraints have a much stronger effect 
when the pt threshold is set at 50 GeV (dashed curve) than at 20 GeV (solid curve); in the 
latter case more phase space is available to radiate gluons. We also show for comparison 
the decorrelation for dijet production at the Tevatron for pj- > 20 GeV. There we see that 
the lower collision energy (1.8 TeV) limits the allowed rapidity difference and substantially 
suppresses the decorrelation at large A. The larger center-of-mass energy compared to 
transverse momentum threshold at the LHC would seem to give it a significant advantage 
as far as observing BFKL effects is concerned. 

The lower set of curves in Fig. ^ refer to Higgs production and will be discussed below. 



3.2 Forward jets from Higgs production via boson fusion 

At the LHC, Higgs bosons can be produced in the process qq qqH, where W or Z 
bosons emitted by the initial state quarks annihilate to form a Higgs [§, see Fig. 1(b). This 
mechanism may be important for discovering a heavy Mh ~ 1 TeV Higgs, since the rate there 



is comparable to the rate from the 'standard' gg ^ H process, see for example Ref. |12]. For 
lighter Higgs bosons, weak boson fusion can provide a complementary method of checking 
that the couplings to heavy quarks and weak bosons are related in the way specified by the 
Standard Model. 

Our interest here, however, is in the final-state quark jets accompanying the Higgs. The 
ability to 'tag' these jets ^ is necessary to suppress backgrounds, for example from qq 
WW, ZZ, and ti production. It is the t-channel W and Z propagators (see Fig. 1(b)) that 
to a large extent determine the regions of phase space populated by the tag jets. Typically, 
they have (9 (TeV) energies and transverse momenta less than M^y. They therefore naturally 
have large (roughly equal and opposite) rapidities. 

Higgs production via WW and ZZ fusion therefore automatically provides a 'BFKL- 
like' dijet sample with large rapidity separation. Figure |^ shows the differential cross section 
(P(y" /dyidy2{yi = -y2 = A/2) as a function of A, for Mh = 500 GeV and pr > 20, 50 GeV 
at the LHC. Notice the difference in shape compared to the QCD LO dijet distributions 
shown in Fig. H the tag jets are naturally produced with a sizeable rapidity separation. 
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The big difference between Higgs and QCD dijet production is of course the t-channel color 
singlet exchange property of the former. This gives rise to a rapidity gap, rather than a 
region between the jets populated by BFKL soft gluons, Fig. l(a).0 The large rapidity 
separation decorrelation between the forward jets induced by multigluon emission between 
them is therefore absent. Interestingly, however, for Mh ^ PrOet) the centrally produced 
Higgs boson acts as a jet decorrelator, i.e. the qq — Hqq matrix element does not depend 
on the relative orientation of pj-i and pt2- Indeed the only significant correlation occurs 
when the tag jets are produced centrally (small A) and all invariants Pi ■ pj are large and 
of the same order. Three-body kinematics then yields approximately back-to-back jets. 
The situation is summarised in Fig. |^, which shows the A dependence of (cosA0) for the 
qq — >• qqH process (lower set of curves). We have again taken Mh = 500 GeV, but in fact the 
dependence on the (heavy) Higgs mass is very weak. Evidently the jets accompanying Higgs 
production are much more decorrelated in azimuth than the Mueller-Navelet QCD jets. In 
principle, therefore, the Higgs tag jets could act as a benchmark for the decorrelation, since 
the theoretical predictions shown in Fig. ^ should be reliable. 



4 Collision energy dependence of dijet production at 
the LHC 

The BFKL increase in the dijet subprocess cross section with A is swamped by the decreasing 
parton distribution functions, as can be seen in Fig. |[ One way 0, H to account for this is 
to take a ratio of cross sections at different A but at the same values of parton momentum 
fraction x so that the pdf dependence cancels out. This requires (at least) two center-of-mass 
energies. The difficulty with this method, as discussed in , is that the cancellation of pdfs 
only occurs using lowest-order kinematics; like any other quantity, the cross section ratio is 
subject to corrections due to kinematic effects. As shown in these kinematic effects turn 
out to be very important at the Tevatron (for -y/i = 630 and 1800 GeV) and the improved 
BFKL MC prediction differs qualitatively from that of naive BFKL. 

In this section we investigate the dijet cross section ratio at the LHC, for the design 
center-of-mass energy ^/s = 14 TeV and for ^/s = 10 TeV, which has been mentioned as a 
possible first step toward the design energy. At the LHC, the higher center-of-mass energy 
compared to transverse momentum threshold may mitigate the strong kinematic suppression 
that spoiled the naive BFKL ratio at the Tevatron. 

We consider a cross section ratio R12 at center-of-mass energies ^/sl and ^/s2 measured 
at fixed values of rapidity difference Ai and A2 chosen so that the asymptotic lowest-order 
QCD cross sections (including pdfs) are equal 0. Given Ai, A2 is given by 

cosh(|Ai) ^ ^ 

COsh(iA2) y/s^ 

''The existence of this rapidity gap has been suggested as an additional signature for Higgs production, 
see for example Ref. | p^ . 
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The cross section ratio to be measured is then 



_ da(v/i„ Ai) 

dcr( VSg' ^2) 

One advantage of this definition is that it uses the rapidity differences between the two 
outside jets, rather than the parton momentum fractions for the event, which can be more 



difficult to measure ||l4| 



In Fig. 1^ we show R\2 at the LHC with ^/s^ = 10 TeV and ^/s^ = 14 TeV as a function 
of Ai, with A2 determined by Eq. (^. Note that although we show results for all values 
of A, it is really only at the larger values that we expect BFKL behavior to be manifest. 
In asymptotic LO QCD Ru = 1 by construction; it is shown for reference as a dotted line 
in the figure. We show the naive BFKL (lower curves) and improved BFKL MC (upper 
curves) ratios for px > 50 GeV (dashed curves) and px > 20 GeV (solid curves). As at the 
Tevatron, we see that the BFKL MC predictions deviate strongly from those of naive BFKL 
due to nonleading effects. The situation is slightly better than at the Tevatron, where the 
cross section was never larger at 1800 GeV than at 630 GeV, in marked contrast to the naive 
BFKL expectation. Again, the LHC benefits from the larger available phase space. 



5 Conclusions 

In summary, the ability of the LHC general purpose detectors to measure 'Mueller-Navelet' 
forward jets with large rapidity and modest transverse momentum will allow an important 
test of QCD 'BFKL' physics. In this paper we have presented the results of an improved 
BFKL Monte Carlo calculation for dijet production at large rapidity separation at the LHC. 
In this approach, subleading effects such as kinematic constraints and the running of the 
strong coupling constant as are incorporated by unfolding the implicit BFKL sum over 
emitted gluons to make the sum explicit .0 

We found that, as at the Tevatron @, 0, these subleading effects can be substantial and in 
particular lead to some suppression of gluon emission. This kinematic suppression is however 
not as dramatic at the LHC, because of its relatively higher center-of-mass energy compared 
to the jet transverse momentum threshold. As a result, improved BFKL MC predictions 
tend to retain more BFKL-type behavior because of the greater phase space for emitting 
gluons. We saw that as the p^ threshold was lowered, the BFKL MC results became more 
naive-BFKL-like. This was true for both the azimuthal decorrelation and the cross section 
ratio at different collision energies. 

Noting that forward jets are also produced at the LHC in heavy Higgs boson production 
via gauge boson fusion, we compared the two dijet production mechanisms. The Higgs case 
is characterized by a lack of hadronic activity in the rapidity region between the two jets 
(the 'rapidity gap'), in contrast to the BFKL case where the region between the jets is 
populated by emitted gluons which give rise to an azimuthal decorrelation of the outside 
jets. Interestingly, however, we found that the Higgs boson itself acts as a decorrelator and 

^Another potentially important effect, not included in the present study, may come from the next-to- 
leading order perturbative corrections to the BFKL equation, which have recently been calculated 
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the jets produced in association with a heavy Higgs show a stronger decorrelation than the 
BFKL jets. 

The LHC holds promise for studies of dijet production in BFKL physics, and such studies 
will benefit from jet detection capabilities that extend far into the forward region and allow 
for pt thresholds as low as possible. However, kinematic effects can still be of quantitative 
importance, and should therefore be incorporated in the theoretical analyses. 
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Figure 1: Schematic representation of dijet production at large rapidity separation A in 
hadron-hadron collisions in (a) the BFKL approach and (b) Higgs production via boson 
fusion. 
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LHC, Vs = 14 TeV 




Figure 2: BFKL and asymptotic QCD leading-order dijet production cross sections at the 
LHC {^/s = 14 TeV) as a function of the dijet rapidity separation. The pdfs are the CTEQ4L 



set |Tll|. The three curves at each transverse momentum threshold use: (i) improved BFKL 
MC with running as (solid lines), (ii) naive BFKL (dashed lines), and (iii) the asymptotic 
(A ^ 1) form of QCD leading order (dotted lines). 
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Figure 3: The azimuthal angle decorrelation in dijet production at the Tevatron {^/s = 
1.8 GeV) and LHC {y/s = 14 TeV) as a function of dijet rapidity difference A. The 
upper curves are computed using the improved BFKL MC with running as', they are: (i) 
Tevatron, pt > 20 GeV (dotted curve), (ii) LHC, pr > 20 GeV (sohd curve), and (iii) 
LHC, Pt > 50 GeV (dashed curve). The lower curves are for dijet production in the process 
qq —>■ qqH for px > 20 GeV (solid curve) and pr > 50 GeV (dashed curve). 
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LHC, qq ^ qqH 




Figure 4: The cross section for the production of jet pairs in association with a Higgs boson 
via gauge boson fusion (Fig. [l|(b)) at the LHC (^i = 14 TeV) as a function of rapidity 
separation of the jet pair for px > 20 GeV (sohd curve), and px > 50 GeV (dashed curve). 
The pdfs are the CTEQ4L set [|lT]. 
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Figure 5: The ratio R12 of the dijet cross sections at the two colhder energies ^/s^ = 
10 TeV and ^/s2 = 14 TeV, as defined in the text for px > 20 GeV (sohd curves), and 
Pt > 50 GeV (dashed curves). The curves are calculated according to: (i) the improved 
BFKL MC predictions using CTEQ4L [|lll pdfs (upper curves), with /i = = 20 GeV, (ii) 
the naive BFKL prediction (lower curves), and (iii) the asymptotic leading-order prediction 
(dotted curve) R12 = 1. 
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